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METHOD AND APPARATUS FOR MULTI-DIMENSIONAL SHAPE 
REPRESENTATION VIA SHOCK FLOWS 



STATEMENT OF GOVERNMENT INTEREST 

1 This invention was partially funded by the Government under grants from the 
National Science Foundation, Contract no. IRI-9700497 and Contract no. IRI-0083231. 
The Government has certain rights in portions of the invention. 

CROSS REFERENCE TO RELATED APPLICATIONS 

2 This application claims the benefit of priority to U.S. Provisional Application 
Serial No. 60/419,563 filed October 18, 2002 and entitled "Three-Dimensional Shape 
Representation Via Shock Flows," the teachings of which are hereby incorporated by 
reference in their entirety. 
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BACKGROUND OF THE INVENTION 

Understanding an object's shape in two or three dimensions is useful in many 
applications including pattern analysis and machine intelligence. For example, 
animations typically employ shape structures upon which forces and feedback are 
exercised. Machine tooling, computer graphics, bio-medical imaging, and architecture 
simulators are further examples where the understanding and efficient representation of 
shape are important. 

Examples of traditional shape models are voxel (that is, volume pixel) or triangle- 
based rendering engines. Other conventional shape models include implicit polynomials, 
generalized cylinders, superquadrics and splines. Another conventional model is the 
Medial Axis (MA), the locus of centers of maximal balls in contact with input data where 
each item of input data is referred to as a generator. Each locus has an associated radius 
value giving the distance from the MA to the generators of that locus. A complementary 
and equivalent definition of the MA is the model that is generated when waves are 
propagated from generators and quenched at loci which trace the MA. In general, the 
MA has a branching structure related to the object's topology and geometry. The MA 
typically encodes a time of formation of each MA point by carrying relative width 
properties via a radius function. The advantages of the MA include that the MA is an 
intuitive representation of elongated objects such as an anthropomorphic form. Further, 
an MA can generally encode the blobbyness of a shape, that is, the varying width of a 
form because the MA has the radius function. Still further, the MA typically makes 
explicit object contour features such as curvature extrema and ridges through the branch 
"tips" of the MA model. The MA describes both the interior and exterior regions of 
space and segments these regions as a function of their relative closeness to the object's 
outline. The MA can partition a shape by combining width and elongation properties, 
where, for example, different MA branches represent different object parts. In general, 
the MA has a hierarchy of scales in that its spatial and time-like properties enable smaller 
features to be distinguished from larger features and the features can be ranked 
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accordingly. Another advantage of the MA is that it provides a fairly complete 
representation. That is, generally, the entire object and its details are represented, from 
small bumps along a boundary to large engulfings of main parts. The MA also employs 
some dimensional reduction, that is, it maps an object and the space that the object 
occupies into a thin set. 

There are a number of conventional methods of creating an MA model of an 
object, also called extracting MA symmetries in 3D. A first conventional method of 
creating an MA is by thinning an object by iteratively peeling off discrete layers of the 
object until an approximation of MA loci remain. A second conventional method of 
creating an MA involves following ridges on distance maps. A third conventional 
method of creating an MA uses "grass-fire" transforms in which "fires" are initiated at 
boundary loci. The fire wavefronts meet and are quenched. The quenching wavefronts 
are the MA loci. A fourth conventional method of creating an MA is to take families of 
primitive shapes which can retro-fit to object data. A fifth conventional method of 
creating an MA is based on using a Voronoi diagram to approximate the MA. A sixth 
conventional method of creating an MA involves doing full bisector computations 
followed by trimming operations to define generalized descriptions of the MA. 

In sum, the MA is a well-developed model of multi-dimensional shapes having a 
variety of methods of making the MA model. 

SUMMARY OF THE INVENTION 

Unfortunately, there are shortcomings to the conventional methods and apparatus 
for representing multi-dimensional objects. Conventional multi-dimensional 
representation technology does not provide a stable, full representation combined with 
fast processing time and data compression. Conventional object representation techniques 
including the MA do not provide robustness against perturbations, partial occlusions and 
articulated movements, choice of metric for matching applications, categorization 
sensitivity, i.e. low-cost matching within a class of similar objects (such as trees) vs. high 
cost matching across classes (such as trees vs. airplanes). MAs for similar objects may 
not be unique and therefore object matching of similar object using MAs is unreliable. 
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Further, the MA is unstable in certain situations for example, for small protrusions. The 
MA typically generates a spurious branch with the introduction of a small protrusion. 
The MA datasets tend to be large which slows down processing time and creates a 
problem for storing MAs. It is desirable to have a method and apparatus for obtaining as 
exact a multi-dimensional representation of an object as possible with a minimum of 
computational limitations. 

Representation of multi-dimensional shapes is accomplished by embodiments of 
the present invention including two computational schemes for a shock scaffold. The 
shock scaffold provides a representation of a multi-dimensional object by linking special 
points of the object called shocks. The shock scaffold is a directed graph that 
incorporates the notion of flow through the links to provide complete information about 
the object shape. The first computational scheme is a Lagrangian scheme and does not 
have a reference grid. The second computational scheme is Eulerian and employs a 
reference grid. Both computational schemes are capable of finding shocks in 
unorganized points clouds, such as sample points generated by laser scanning of the 
object. Both computational schemes involve generating wavefronts from shocks to 
determine the organization of the shape of the object. 

The first computational scheme for obtaining a shock scaffold is the Lagrangian 
computational scheme which does not require a fixed grid. In this method, in initial set 
of sources for the flow along the MA is identified. These are sources for MA sheets, that 
is, pairing two distinct outline pieces. The sources are then paired to identify sources of 
flow for MA curves. The MA curves are then paired to identify MA vertices where the 
curves implicitly intersect. Some of the sources of flow for curves and some of the 
vertices are relays for the flow, that is, loci where the entrant flows collide and are re- 
initiated, in different directions. The new flows are tracked and the system looks for new 
intercepts (pairings) until all flows terminate at sinks, or at infinity, or at the spatial limits 
of a fixed zone where all computations are restricted to take place. 

A second computational scheme for obtaining a shock scaffold is the Eulerian 
grid. Shocks are discovered by the collision of waves propagated in a tessellated space 
which takes the form of a rigid uniformly connected grid. An explicit search is 



Attorney Docket No.: BRU02-01 



-5- 

performed through the embedding space. Essentially, the propagation carries only those 
samples to the grid node where a pairing must occur, thus explicitly limiting the number 
of false pairings. The computation complexity of the Eulerian method is essentially 
dependent on the number of grid nodes, or cells, used to tessellate space. 

1 1 The present invention operates to provide an object model even where the initial 
data is a collection of unorganized points or polygonal patches in space, e.g., unorganized 
point clouds. These are arbitrary sets of points in three-dimensional space for which no 
initial connectivity information is known. Ability to operate using only point clouds as 
input has applications because objects conventionally are sampled by points. Many data 
collection sensors, such as, active sensing by laser scanners, produce point samplings. 
Another advantage in capability of analyzing points clouds is that connectivity data, 
where it is provided, can be unreliable or incomplete. The ability to operate in an 
environment where some links between points are incorrect and where some needed 
links are missing makes the system of the present invention more reliable and functional 
than conventional systems. 

12 The system is also capable of modeling objects where the input data is surface 
patches also rather than simply points. The system according to the present invention is 
therefore able to use data stored as surface meshes. The surface patches can be 
unorganized polygonal clouds. 

13 More specifically, embodiments of the invention provide methods and apparatus 
that receives as input unorganized points in space that are sampled from the surface of an 
object and derives a shock scaffold representation of the object. The method forms a 
plurality of shocks based on the unorganized point input where the shocks are formed at 
collisions of a plurality of wavefronts initiated from the plurality of shocks. The shocks 
hold topology information about the surface of the object including flow speed and 
direction from the boundary of the object. The method then generates the whole shock 
scaffold from the plurality of shocks. The object's surface can be reconstructed from 
data in the shock scaffold representation. 

14 In another embodiment of the invention, the step of forming shocks from the input 
further comprises the step of defining a plurality of clusters where the input data is 
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distributed among the defined clusters. The method then examines each cluster to 
determine pairs of generators by applying visibility constraints to the generators. A 
shock candidate is generated from each pair of generators. The method then validates the 
shock candidates by examining the contact sphere of each candidate. If no other 
generators other than the generator pair are in the contact sphere, then the shock 
candidate is a valid shock. If there are other generators in the contact sphere, then the 
shock candidate is not valid. In this step, shocks defining shock sheets are computed. 
These types of shocks are used to find shock curves which are in turn used to find shock 
vertices of the shock scaffold. In another embodiment of the invention, clusters of 
generators are treated as meta-generators. At the cluster scale, a first cluster acts on a 
second cluster similarly to the way a first generator acts on a second generator in a 
generator pair. The examination of cluster pairs enables the system of the invention to 
rapidly rule out sections of space from being searched for other cluster pairs. At the 
generator scale, the content of a cluster is considered by pairing generators within a 
cluster. This method of fine-scale generator pair selection may be effectively used in 
dynamic data acquisition systems because the method computes constrains for points in 
clusters as the points are acquired, as though from a dynamic data acquisition system. 

In another embodiment of the invention, the step of computing shock candidates 
further comprises defining a fixed multi-dimensional grid in space around the plurality of 
points. The fixed grid acts as a reference for wavefront propagation from each of the 
plurality of points. The method initiates logical wavefronts by initiating cellular 
automata along the grid from chambers including at least one of the plurality of points. 
The cellular automata collide. The method detects a shock at each collision. In this 
alternative method of finding shock candidates, the fixed grid provides a reference and an 
initial approximation of organization of the input data. The use of a fixed grid permits 
the method to avoid validating shocks in an exhaustive manner in order to ensure the 
associated ball is maximal (i.e. empty of other generators). 
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BRIEF DESCRIPTION OF THE DRAWINGS 

16 The foregoing and other objects, features and advantages of the invention will be 
apparent from the following description of particular embodiments of the invention, as 
illustrated in the accompanying drawings in which like reference characters refer to the 
same parts throughout the different views. 

17 Figure 1 is a block diagram of a computer system including a shock scaffold 
engine according to principles of the invention; 

1 8 Figure 2 is a sketch of a plurality of sampled points of a box-shaped object; 

1 9 Figure 3 is a shock scaffold generated from the sampled points of Figure 2; 

20 Figure 4 is a diagram of a plurality of objects and the medial axis representations 
of each of the plurality objects and shock structures of each of the plurality of objects; 

21 Figure 5 is a diagram of a two-dimensional object illustrating shock generation; 

22 Figure 6 A is a medial axis representation of a three-dimensional object; 

23 Figure 6B is a shock hypergraph representation of a three-dimensional object; 

24 Figure 6C is an augmented shock scaffold representation of a three-dimensional 
object; 

25 Figure 6D is a shock scaffold representation of a three-dimensional object; 

26 Figure 6E is a reduced shock scaffold representation of a three-dimensional 
object; 

27 Figure 6F is a topological scaffold representation of a three-dimensional object; 

28 Figure 7 is a side view diagram of contact spheres intersecting boundaries giving 
rise to types of shocks; 

29 Figure 8 A is a diagram of an overview of shock scaffold computation using a first 
computational method according to principles of the present invention; 

30 Figure 8B is a diagram of an overview of shock scaffold computation using a 
second computational method according to principles of the invention; 

3 1 Figure 9 is a flow chart of the general operation of the first computational method 
for deriving a shock scaffold according to principles of the invention; 

32 Figure 10 is a flow chart of the operation of shock candidate computation of the 
first computational method of Figure 9; 
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33 Figure 1 1 is a flow chart of the operation of validating shock candidates computed 
using the method of Figure 10; 

34 Figure 12 is a diagram illustrating visibility of generators according to principles 
of the invention; 

35 Figure 13 is a diagram illustrating paired generators and the associated dead zone 
and visible zone of the paired generators according to principles of the invention; 

36 Figure 14A is a diagram illustrating a one-dimensional beam according to 
principles of the invention; 

37 Figure 14B is a diagram illustrating a two-dimensional beam according to 
principles of the invention; 

38 Figure 14C is a diagram illustrating a three-dimensional beam according to 
principles of the invention; 

39 Figure 1 5 is a two dimensional grid showing beam propagation using the Eulerian 
computational method of deriving a shock scaffold of the present invention; 

40 Figure 16 is a flow chart of the basic operation of the Eulerian computational 
method for deriving a shock scaffold of the present invention; 

41 Figure 1 7 A is a diagram of a grid chamber 920 illustrating passage of the first 
type of Lagrangian cellular automata, the 2D Lagrangian automata for shock sheets 
according to principles of the invention; 

42 Figure 17B is a diagram of a grid chamber 950 illustrating passage of the second 
type of Lagrangian cellular automata, the ID Lagrangian automata for shock curves 
according to principles of the invention; 

43 Figure 18 is flow chart of the operation of the Lagrangian automata in finding 
shock sheets and shock curves in deriving a shock scaffold according to principles of the 
invention; 

44 Figure 19A is an example of a Voronoi diagram of a set of point generators which 
is suitable for use with the computational methods of the present invention; 

45 Figure 19B is a medial axis for the set of point generators of Figure 19A; and 

46 Figure 19C is the shock scaffold for the point generators of Figure 19A. 



Attorney Docket No.: BRU02-01 



-9- 

DETAILED DESCRIPTION 
Introduction 

47 Representation of two and three-dimensional shapes is accomplished by 
embodiments of the present invention including two computational schemes for a shock 
scaffold. The shock scaffold provides a representation of a multidimensional object by 
linking special points of the object called shocks. The shock scaffold is a directed graph 
that incorporates the notion of flow through the links to provide complete information 
about the object shape. Flow is the projection of the gradient of the radius function 
associated with the medial axis (MA) maximal contact balls. This projection is taken 
onto the tangent space to MA sheets or curves. The projection creates a vector field 
along the MA sheets and curves. Certain singularities of this induced vector field are 
then retained to constitute the nodes of the shock scaffold. The connectivity of the 
singularities via the flow creates the directed links of the scaffold. In short, the shock 
scaffold enables information from an N-dimensional space to be mapped to a graph while 
the MA maps the information to a set of N-l -dimensional sets, e.g. surfaces in three- 
dimensional space. 

48 The first computational scheme is a Lagrangian scheme and does not have a 
reference grid. The second computational scheme is Eulerian and employs a reference 
grid. Both computational schemes are capable of finding shocks in unorganized points 
clouds, such as sample points generated by laser scanning of the object. Both 
computational schemes involve generating wavefronts from shocks to determine the 
organization of the shape of the object. 

49 The surface of an object is extracted from a shock scaffold by ranking selected 
MA curves on the basis of their associated local surface interpolant. Each MA curve is 
generated from three associated input points, and the local surface interpolant associated 
with the MA curve is defined as the triangle interpolating the three points. Once the 
surface is recovered, the shock scaffold is segregated into two subgraphs, one linked to 
the sampled surface and other lying on either side of the recovered surface. The non- 
surface sub-graph is connected in larger surface and curve structures resulting in an 
approximation of large MA sheets and curves lying deeper inside and out of the 
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originally sampled object. In the case where the goal is surface meshing alone, the 
calculation of the scaffold may be ended at the surface interpolant. Once the MA curves 
providing the surface mesh are found, the remaining shock nodes need not be calculated 
if the goal of calculation was only the surface. 
Representation of Shape 

50 Figure 1 is a block diagram of a computer system suitable for use by the present 
invention. The computer system 100 includes a controller 105 having a shock scaffold 
engine 1 10, a memory 1 15, a data store 120, and a data interface 125. The computer 
system 100 further includes a display 130 and a data collection system 135. 

5 1 The shock scaffold engine 1 10 is part of the controller 105 and acts in conjunction 
with the memory 115 and data store 120 in the computer system 1 10 to compute shock 
scaffolds from input data about an object such as sampled points from a laser scanner. In 
this embodiment, the sample points are received from the data collection system 135 
connected to the computer system 100 through the data interface 125. The data 
collection system 135 may, for example, be a laser scanning system to scan object 
surfaces. The shock scaffold is a data structure stored in the data store 120 and can also 
be visually displayed on the display system 130 connected to the computer system 100. 

52 Figure 2 shows a set of sampled points 155 which is suitable for use by the 
computer system 100. By way of example only, the sampled points 155 are points from a 
surface of a box-shaped object 150, and sampled points 155 from other shapes are 
suitable for use by the computer system 100 as well. As will be explained in further 
detail shortly, the shock scaffold engine 1 10 of the computer system 100 is capable of 
obtaining the sampled points 155 through the data collection system 135, generating a 
shock scaffold representing the box-shaped object 150 from the sampled points 155, and 
storing the shock scaffold in the data store 120 in an efficient amount of memory space 
for subsequent use (e.g., for rendering images of the box-shaped object 150 on the 
display 130 based on the shock scaffold when the shock scaffold is stored in the data 
store 120, also see Figure 1). 

53 Even though the sampled points 1 55 may be recognizable to the human eye as 
being from the surface of the box-shaped object 150, it should be understood that there is 
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no linking information included with the sampled points 115. Accordingly, conventional 
computers generally view the sampled points 155 as an unorganized point cloud, and are 
unable to classify the sampled points 155 as being from the box-shaped object 150. In 
contrast to such conventional computers, the computer system 100 is capable of 
analyzing the sampled points 155 to recognize the surface of the box-shaped object 150, 
as well as performing similar operations on more complex data. As a result, the 
computer system 100 is well-suited for a variety of applications in which a set of sampled 
points but no linking information is available (e.g., medical imaging, archeology, etc.). 
No assumption is made on the input points (which may be degenerate configurations) and 
the sampling methods (which need not be uniform)." While the box-shaped object 150 
shown in Figure 3 is a simple shape, more complex objects such as a human brain or a 
pot shard, can be represented multi-dimensionally including the object's edges and 
surface details. 

54 Figure 3 shows a shock scaffold 200 generated by the computer system 100 from 
the sample points 155 of Figure 2. In general, each shock scaffold 200 generated by the 
computer system 100 includes a set of shock curves 205 and a set of shock nodes 210. 
The location of the shock nodes 210 within the shock scaffold 200 is with respect to the 
object's boundary, in this case, the boundary of the box-shaped object 150. A shock node 
210 is a point augmented with attributes about the flow. A shock curve 205 is a link 
between shocks. A sheet 215 is defined by a plurality of shock curves 205 defining a 
loop implicitly defining a boundary for the sheet. The sheet 2 1 5 therefore also connects a 
plurality of shock nodes 210. Together, the shock nodes 210, shock curves 205 and 
sheets 215 provide a full representation of the box-shaped object 150. The flow in a 
shock scaffold is induced by the gradient of the radius function of the shock scaffold. 

55 It should be understood that the shock scaffold 200 is a data structure capable of 
being represented in the computer memory 1 15 and stored as data in the data store 120 as 
well as a visual representation of a multi-dimensional object capable of being displayed 
on a visual display 130. The shock nodes 210, shock curves 205 and sheets 215 are 
described herein as physical entities but it should be understood that the shock nodes 210, 
shock curves 205 and sheets 215 are signals generated in the computer system 100. 
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While the shocks nodes 210 are the main data upon which the scaffold is built, the 
operation of the invention will be explained in terms of how their representations interact 
in order to facilitate understanding of the invention. 

56 Figure 4 is presented to illustrate differences between the MA and the shock 
scaffold and also to illustrate advantages of representing multi-dimensional objects using 
the shock scaffold rather than the MA. Figure 4 shows a plurality of objects 250-1, 250- 
1', 250-2, 250-2', 250-3, 250-3', 250-4, 250-4', 250-5, 250-5' (collectively numbered 250). 
Also shown are the medial axis representations 255-1, 255-2, 255-3, 255-4, 255-5 
(collectively numbered 255) of the objects 250 and the shock graph representations 260- 
1, 260-2, 260-3, 260-4, 260-5 (collectively numbered 260) of the objects 250. A shock 
graph is a shock scaffold in two dimensions. The objects 250 and representations 255, 
260 in Figure 4 are presented to illustrate the concept of flow in a shock structure 
representation of an object and further to illustrate the distinction of shock scaffold from 
the medial axis. 

57 The objects 250-1', 250-2', 250-3', 250-4', and 250-5' contain their shock scaffold 
representations 260. In contrast to the MA representations, each object 250 has an 
essentially different shock graph representation 260 having two or more shock node 210 
connected by curves 285. The shock graph representations 260 incorporate the notion of 
flow indicated by arrows 290, 292. The arrows 290, 292 indicate the direction of the 
flow at the shocks 280 and through the curves 285. Flow has both direction and speed. 
The direction of the shock node 210 is the direction of the flow at the point location of 
the shock node 210. Flow is directed from, through and to shocks in the direction of 
increasing radii values from the boundaries of the object. Speed is indicated in the shock 
structure representations 260 by the number of arrows 290. A single arrow 290 indicates 
slower speed and double arrows 292 indicate faster speed and possibly other attributes, 
such as accelerations. This illustrates an advantage of shock scaffold representation. The 
shock structure representations 260 provide a model of the objects 250 that distinguishes 
one object from another through the incorporation of flow. 

58 Figure 5 is presented to illustrate the relationship between shock nodes 210 
(Figure 3) and flow (represented by arrows) and surface boundaries of an object. The 
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object 300 shown in Figure 5 is a two-dimensional object. Two dimensions are used for 
simplicity and clarity. In general, shock nodes 210 and flow have equivalent 
relationships to object boundaries in three dimensions. Figure 5 also illustrates how 
shock nodes 210 are found from surface boundaries in order to better understand how the 
shock scaffold engine 1 10, using the methods of the present invention, computes shock 
nodes 210 from sampled points of the object's surface. 

59 When an object's surface is already known, the computer system 100 can readily 
generate a shock scaffold 200 representing the object's surface. Any place of an 
extremum of shape has a point that represents the center of the curve along the boundary 
at that place. One way of deriving the curve center 385, which is one type of shock, is by 
finding the center of an osculating disc 340 at the boundary curve at a place of the 
extremum 345. Depending on the orientation of the boundary curve, the osculating disc 
340 may touch the boundary either on the inside or the outside of the object. The 
osculating disc 340 is also called a contact disc. The three-dimensional equivalent of the 
osculating disc is the osculating sphere, also called the contact sphere. A shock structure 
302 arises from wave propagation from the boundaries of the object 300. As before flow 
is induced by the gradient of the radius function associated with the contact discs. 

60 The two-dimensional object 300 has a plurality of points 305, 310, 315, 320, 
325, 330, 335 inside the object 300. Also, inside the two-dimensional object 300 are a 
plurality of graph segments 350, 355, 360, 365, 370, 375, 380 connecting some of the 
points 305, 315, 320, 325, 330, 335 to form a shock graph 302. The points 305, 310, 315, 
320, 325, 330, 335 illustrate the various points within the object that are either shock 
locations or that do not form shocks. The points 305, 310, 315, 320, 325, 330, 335 are 
provided for illustration purposes and are not considered to be comprehensive as there are 
both more points that are shock locations and points that are locations where no shocks 
are formed. There are no shocks at point 310 where no contact sphere can be maximal. 
There is an initial shock at point 3 1 5, a point corresponding to a curvature extremum of 
the boundary. There is a shock at point 305, a junction of two in-coming flows which is 
also called a shock sink. There is a shock at point 320 where there are two outgoing 
flows also called a shock source. There is a shock at point 325 where the flow is 
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monotonic, or regular, which also defines a shock segment, also called a shock relay. 
There is a shock at point 330 where a pair of branches terminate. There is another shock 
at 335 where 3 branches terminate. The shock sinks, sources and relays each provide 
topology information about the two-dimensional object 300. In three dimensions, the 
shock scaffold has flow along sheets and curves and through the vertices of the shock 
scaffold as seen below in Figure 6. The shock structure 302 is formed by connecting the 
shocks 305, 315, 320, 325, 330, 335 and the graph segments 350, 355, 360, 365, 370, 
375, 380 which, here, are shock curves. Where the surface of the object is not known, 
and the input is an unorganized point cloud, the computer system 100 needs to use other 
methods to derive a shock scaffold, for example the methods of the present invention. 

61 Figure 6 illustrates another advantage of the shock scaffold over the MA as a 
representation of multi-dimensional objects beyond the advantages discussed with regard 
to Figure 4. Figure 6 also further defines and describes the shock scaffold with respect to 
the MA. The shock scaffold provides a complete representation of a multi-dimensional 
object while providing greater data compression than the MA. 

62 Figure 6 shows several semi local representation models for representing multi- 
dimensional objects. Figure 6A shows a medial axis representation 400, or MA; Figure 
6B shows a shock hypergraph 405; Figure 6C shows an augmented shock scaffold 410; 
Figure 6D shows a shock scaffold 415; Figure 6E shows a reduced shock scaffold; and 
Figure 6F shows a topological scaffold. 

63 The medial axis representation 400, or MA, of Figure 6A is formed by three 
medial sheets 420 intersecting into a medial curve 425. While the MA 400 is a full 
representation of a multi-dimensional shape, the MA 400 uses a larger amount of data to 
make the full representation than does the shock scaffold 410 because the medial curves 
and medial sheets, as well as the nodes of the MA are explicitly stored. 

64 As shown in Figure 6B, the shock hypergraph 405 is an intermediate 
representation between the MA 400 and a shock scaffold 410. The shock hypergraph 405 
is derived from contact spheres and flow as described above in two dimensions with 
respect to Figure 5. In a shock hypergraph 405, shock nodes 432 are identified and 
connected by directed links 430. The points represented by triangles are shock nodes for 
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shock curves, e.g. a shock source for a curve, as in an A\ -2 . The small squares 
represent shock nodes for shock sheets (e.g. a shock source for a sheet, as in an j£ x -2). 

The shock hypergraph is a complete representation of a shape, but the shock hypergraph 
has some redundancy. The shock hypergraph makes explicit the shock nodes and their 
connectivity in the interior of sheets and curves. 

65 As shown in Figure 6C, the augmented shock scaffold 408 is a representation 
where shock nodes along curves are connected by directed links. The triangles represent 
shock nodes for shock curves. The cycle order of the hyperlinks is indicated by a 
counterclockwise arrow in each sheet. The augmented shock scaffold explicitly traces 
the interiors of shock sheets, but without shock nodes for these sheets. 

66 As shown in Figure 6D, the shock scaffold 410 is a shock hypergraph without the 
redundancy. The shock scaffold 410 has three intersecting shock sheets 215. Each shock 
sheet 215 is defined by shock nodes 210, some of the shock nodes 210 forming 
boundaries 445 to the sheets 215. The circles represent shock vertices. The triangles 
represent shock nodes for shock curves. The interiors of the shock sheets 215 are 
implicitly defined by the shocks 210 and boundaries 445. Another definition of a shock 
scaffold, therefore, is a hierarchical data organization where either shock sheets 215 or 
shock sheets 215 and curves are represented by their boundaries 445. Because the shock 
sheets 215 are implicit, the shock scaffold 410 provides an object representation with 
greater data compression than the MA 400. 

67 As shown in Figure 6E, in the reduced shock scaffold 415, the shock sheets 215 
and boundaries 445 (which are also shock curves) are implicitly defined by the shock 
nodes 210. A reduced shock scaffold 415 provides an object representation using a 
minimal set of data because both shock sheets 215 and shock curves are implicit. In 
some instances, the reduced shock scaffold 415 does not provide a complete 
representation of an object. 

68 As shown in Figure 6F, in the topological scaffold 418, only the topology of the 
graph structure of the shock scaffold is retained. The topological scaffold 418 is a set of 
nodes and links without geometry. The topological scaffold is useful for defining larger 
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or coarser classes of objects and for applications where only the connectivity (i.e., the 
topology) is important. 

69 There are different types of shock nodes 210 and it is important to understand 
some of the types of shocks as the computational methods presented herein generate a 
shock scaffold 200 by identifying shock structures within unorganized input data. As 
discussed above with regard to Figure 5, the center of a maximal disc (2D) or maximal 
sphere (3D) is a shock nodes 210. The shock nodes 210 has a type according to the type 
of contact of the maximal sphere with the object surface and its flow type. The 
computational methods of the present invention generate a shock scaffold 200 from input 
data about an object's surface by identifying shock structures in the input data as defined 
by particular shocks 

70 Figure 7 is presented to describe some of the types of shock nodes 210 and also to 
present nomenclature for the various types of shocks. Figure 7 is a side view of contact 
spheres 480 giving rise to some types of shocks 482, 484, 486, 488 and 490. Shocks are 
first classified according to the number of contact boundary points, and according to the 
degree of contact A shock J^ k has a contact sphere touching a boundary at n points, 
each point having a k+1 degree of contact. 

71 Table 1 shows a classification of the 1 8 possible shock points based on contact 
with spheres, jf k , and flow type. The later flow typology is explained next. 
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72 First flow typology along sheets and curves is considered. A regular shock, also 
referred to as a 1 st order shock is a shock at which flow goes through smoothly: (i) along 

a sheet: (") along a curve: J*\ ~ ^ A3 ~ 1 ■ A shock source > also re f erred t0 as a 

2 nd order shock is a shock which initiates flow: (i) on a sheet: Al\ - 2 ; (ii) on a curve: 

j£ x - 2, A?, ~ 2 • A shock relay, also referred to as a 3 rd order shock is a shock which is 

both a source and a sink for the flow: (i) for a sheet: j[ x - 3 ; (ii) for a curve: 

- 3, - 3 . A shock sink, also referred to as a 4 th order shock is a shock at which 

flow type terminates: (i) for a sheet: ^ - 4 ; (ii) for a curve: A\ ~ 4 ' A3 ~ 4 • Second > 
the typology of flow at vertices is based on the number of inward flows, i.e., along shock 
curves intersecting at the vertex. 

73 Returning to Figure 7, for shock 482, the contact sphere 480 touches a boundary 
500 of an object at one contact point 502. The degree of contact k = 1, so shock 482 is 
classified as an^ type of shock. Where n=l, the number of contact points not shown by 
convention. For shock 484, the contact sphere 480 touches the boundary at one contact 
point 504 with a degree of contact k=2. Accordingly, shock 484 is classified as an A 2 

type of shock. For shock 486, the contact sphere 480 touches the boundary 500 at one 
contact point 506 with a degree of contact k=3. Accordingly, shock 486 is classified as 
an A 2 tyP e of shock. For shock 488, the contact sphere 480 touches the boundary 500 at 
two contact points 508, 510 with a degree of contact k=l . Accordingly shock 488 is 
classified as an j£ type of shock. Sets of £ shock define the interior of sheets within 

the shock scaffold Finally, for shock 490, the contact sphere 480 touches the boundary 
500 at three contact points 512,514,516 with a degree of contact k= 1 . Accordingly, 

shock 490 is classified as an Al tyP e of snock - 

74 Only odd orders of contacts correspond to maximal contact spheres (i.e. where the 
surface outline does not cross over the boundary of the maximal contact sphere). Once 
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the contact typology is defined, the associated radius flow typology is determined. That 
is, the gradient of the radii of the spheres is projected in the tangent space of the medial 
sheets and curves. The shock scaffold nodes are identified from the singularities of the 
induced vector field. This yields the sources for sheets and curves, relays for sheets, 
curves and at vertices, and sinks for sheets, curves and at vertices. All vertices, being 
points, are used as nodes of the scaffold. Three types of vertices are distinguished from 
the number of outward flows at the particular vertex. Referring to Table 1 , the first 
column (regular shock flow) is made implicit in the scaffold. The corresponding MA 
points represent the interior of sheets and curves where no singularities of the flow occur. 
Degenerate cases, where an entire MA sheet patch (curve or surface), or an entire MA 
curve segment is produced at once, are labeled as special types of relays (numbered "3" 
in analogy to other relays) and represented by a single arbitrary point along their trace. 
The singularities of the flow being isolated nodes in space are then connected by 
following the flow direction to create the directed graph structure that is the basis of the 
shock scaffold hierarchy. 

Where the generators are points or spheres, many shock types cannot occur. This 
simplifies the typology of shocks to eight shocks of the sixteen shocks of Table 1 . These 

eight possible shock points based on contact with spheres, j£ , and flow type, when 

generators are points only, are shown in Table 2. When studying only points and spheres, 
those shocks occurring only due to ridges can be eliminated from consideration. 
Accordingly, those shocks in the "rib" and "rib end" rows of Table 1 are eliminated from 
consideration in this case. Further, the shock scaffold in this case includes the Voronoi 
diagram as a special case. The system of the present invention therefore permits the 
efficient computation of a three dimensional Voronoi diagram as a by-product, as will be 
further explained below. 
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Computational Aspects 

76 Figures 8A and 8B respectively illustrate generally the Lagrangian and the 
Eulerian methods of shock scaffold computation of the present invention. Figure 8A is a 
diagram of Lagrangian computation. Points Gi, G2, and G3 are three points of data 
sampled from the boundary of an object. The method of tracking centers of contact 
spheres described above with regard to Figure 5 is computationally intensive, and also 
not generally practical where the boundary of the object is not already known. Instead, in 
the present invention, the shock scaffold engine initiates a logical wave at input elements, 
also referred to as generators, such as points Gj, G 2 , and G 3 along the boundary of the 
object being represented. The waves meet at shock points 550. The shock points are 
shocks that define shock sheets within the shock scaffold of the object from which the 
generator points were sampled. In the Lagrangian method, the input elements are paired 
to create initial shock sources defining sheets representative in the shock scaffold being 
derived. These shock sources of sheets are then paired to create shock sources for curves. 
An iterative process of pairing detected sheets and curve shock representatives then 
permits the construction of the remaining scaffold shock nodes and the connectivity 
between these. 

77 Figure 8B is a diagram of the Eulerian computation. In the Eulerian method, 
input elements such as points Gi, G2, and G 3 are located with respect to a fixed grid 555. 
A propagation from each input element, or generator, is then performed on the fixed grid 
555. The grid nodes closest to each input element are then labeled accordingly. 
Elementary cells of the grid where two or more labels arrive at the "same time" are used 
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to compute shocks. The trace of the shock sheets and curves is traced explicitly by 
computing intercepts with the fixed grid. 

78 In alternative embodiments of either computation method, the input elements are 
surface patches such as triangles or polygons. Additionally, different types of scanning 
processes produce different degrees of object surface information. The data having the 
least amount of information and therefore the greatest degree of difficulty in computing a 
shock scaffold is the unorganized point cloud. The main cause of computational 
difficulty is that without a priori knowledge of connectivity between points in the 
unorganized point cloud, any two points (generators) is a likely pair to create a shock. 
The computational methods described herein are applied to unorganized point clouds, 
however, it should be understood that even greater efficiency of shock scaffold derivation 
can be obtained where more information about data organization is obtained, for example 
points on the object boundary together with the surface normals at those points. 

79 Figure 9 is a flow chart of a procedure performed by the shock scaffold engine 

1 10 (Figure 1) when operating in accordance with the Lagrangian computational method 
for deriving a shock scaffold from an unorganized set of data such as an unorganized 
point cloud. Given a set of points located in three-dimensional Euclidean space, the first 
computational method recovers a shock scaffold efficiently by identifying a local 
network of connectivity between points thereby identifying sources of flow and then by 
propagating flow to recover the remainder of the shock scaffold. 

80 At step 600, the shock scaffold engine 110 receives as input unorganized input 
elements such as an unorganized point cloud. In an alternative embodiment of the 
invention, the unorganized input elements are surface patches of an object. 

8 1 At step 602, the shock scaffold engine 1 1 0 computes shock candidates. In this 
step, the shock scaffold engine 110 examines the input elements to discover sources of 
flow, referred to as shock candidates. The shock candidates provide initial points of 
organization of the shock scaffold. 

82 At step 604, the shock scaffold engine 1 10 validates the shock candidates as will 
be described below. The validated shock candidates are shocks that define shock sheets 
of the shock scaffold. 
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83 At step 606, the shock scaffold engine 1 10 intersects shock sheets to find shock 
curves. The intersection of a sheet with other sheets that share a common generator 
creates shock curves. The shock curves are added to the structure of the shock scaffold 
being derived. 

84 At step 608, the shock scaffold engine 110 uses the shock curves found in step 
606 to find shocks vertices and also new shock sheets. The intersection of a curve with 
other curves sharing two common generators creates new shock vertex candidates. The 
shock vertices and shock curves are added to the shock scaffold structure being derived. 
The shock scaffold engine 1 10 returns to step 606 to find shock curves with the new 
shock sheets. 

85 At step 610, the shock scaffold engine 110 uses the shock vertices found in step 
608 to find new shock sheets and shock curves. The shock scaffold engine 1 10 returns to 
step 606 to find further shock curves with the new shock sheets. The shock scaffold 
engine 1 10 returns to step 608 to find further shock vertices with the new shock curves. 

86 In an alternative embodiment of the invention, the unorganized input elements at 
step 600 is an incomplete set of data about the object to be represented by the shock 
scaffold. An incomplete data set is received as input when, for example, analysis of data 
from a scanner is performed while scanning continues. At step 612, in the alternative 
embodiment of the invention, the shock scaffold engine 1 10 looks for new input 
elements. If new input elements are received, the shock scaffold engine 1 10 returns to 
step 602 to compute new shock candidates from the new input elements. If no new input 
elements are received the shock scaffold engine 1 10 proceeds to end the process. This 
alternative embodiment can be used in dynamic data acquisition scenarios, as described 
above, and in multi-resolution scenarios, where a subset of the data is first processed, and 
known remaining data is successively added an processed toward the finest resolution. 
This method takes advantage of the visibility constraints (described below with respect to 
Figure 10) to optimize the search for valid pairs. 

87 At step 614, the process ends when the shock scaffold engine 110 has iterated 
through all of the shock vertices and there are no more new shock sheets and curves that 
is, when there are no more shock with unlinked outward flows. 
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88 Figure 10 is a flow chart of a procedure to compute shock candidates performed 
by the shock scaffold engine 1 10 (Figure 1) when operating in accordance with the 
Lagrangian computational method for deriving a shock scaffold from unorganized input 
elements, step 602 above. The shock scaffold engine 110 identifies from the unorganized 
point cloud shocks that are candidates for defining shock sheets. 

89 At step 620, the shock scaffold engine 1 10 defines clusters of input elements, 
referred to as generators. In order to efficiently recover a shock scaffold from the 
generators, the shock scaffold engine 110 imposes a coarse structure on the generators so 
that the shock scaffold engine 1 10 has smaller sets of generators to work on and also to 
better discover local relationships among the generators. There are many ways of 
imposing a coarse structure on the generators, such as adaptive bucketing. The input 
elements are defined within Euclidean space and thus are defined in three dimensions 
having an x direction, a y direction and a z direction. In forming clusters using adaptive 
bucketing, the shock scaffold first divides the space containing the generators into slabs 
along the x direction, called x-buckets, where the x-buckets contain approximately equal 
numbers of generators. In some data samplings, some of the x-buckets will be empty. 
The shock scaffold 110 then divides the space in each non-empty x-bucket along the y 
direction to form xy-buckets. The shock scaffold engine 110 forms the xy-buckets such 
that those xy-buckets containing generators contain approximately equal numbers of 
generators. Finally, the shock scaffold engine 1 10 divides the space of each xy-bucket 
along the z direction into xyz-buckets where the xyz-buckets contain approximately equal 
numbers of generators. 

90 Returning to Figure 10, at step 622, the shock scaffold engine 1 10 pairs 
generators within each cluster. Two generators form a pair of generators if they are 

visible to each other. Visible means that the two generators could form an ^ shock 
source. Generators forming pairs are typically close to each other, that is the generators 
are neighbors. Figure 12 is a diagram of visible neighbors. Figure 12 shows three 
generators, G i5 G j? and G k , a first shock 700 and a second shock 702. Generators Gj and 
Gj are contained in a first maximal sphere 704 and generators Gj and G k are contained in 
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a second maximal sphere 706. Gj is visible from G i5 but G k is not. Gj and G k are both 
visible from Gj. Only Gj is visible from G k . Therefore, Gj and G k cannot be paired to 
form a maximal sphere defining a shock sheet because the presence of G k does not allow 

the formation of a maximal sphere with an ^ shock source. That is, Gj is M in the way." 

91 Continuing with step 622, the shock scaffold engine 110 selects a generator 
within a cluster. The shock scaffold engine 1 10 then tests the selected generator with 
respect to the other generators within that cluster using visibility constraints. The 
visibility constraints include finding a closest neighbor to the selected generator. For any 
generator, its closest neighbor is always a visible neighbor. When the shock scaffold 
engine 1 10 finds a generator forming a pair with the selected generator, the shock 
scaffold engine defines a dead zone. The dead zone is the space behind a tangent plane 
of the contact sphere containing the selected generator and the paired generator, tangent 
at the paired generator. 

92 Figure 13 is a diagram illustrating paired generators and their visible zone 720 
and dead zone 722. Figure 13 shows generators G i? Gj, and G k of Figure 12. In this 
example, generators Gj and Gj are paired. The dead zone 722 is the space behind the 
tangent plane 724 of Gj, which is "invisible" to Gj and therefore contains no generators 
with which Gj can be paired. The zone visible to the generator Gi is the space in front of 
the tangent plane 724. Defining a dead zone quickly eliminates generators from 
consideration. 

93 When the shock scaffold engine 1 10 selects a next generator to be tested, the 
shock scaffold engine 110 first determines whether the generator is located within the 
dead zone 722 before applying visibility constraints. A test for whether a generator Gk is 
in the dead zone 722 of generator Gj is by verifying that the following relationship is 
valid: 

94 As generators that pair with the selected generator are found, the dead zone 722 
enlarges and the visible zone 720 decreases because dead zones defined by the shock 
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scaffold engine 1 10 after pairing new generators are additive. In an alternative 
embodiment of the invention, dead zones resulting from invalid pairings of generators are 
retained in order to quickly restrict the visible area around the selected generator. 

95 At step 624, after the elements within each cluster are paired, the shock scaffold 
engine 110 selects a cluster from the defined clusters. 

96 At step 626, the shock scaffold engine 1 1 0 pairs the selected cluster with 
neighboring clusters in a process similar to pairing generators within a cluster. The shock 
scaffold engine 110 applies visibility constraints and defines dead zone using the clusters 
as meta-generators. In an alternative embodiment of the invention, each cluster is 
represented by one or more virtual generators to reduce complexity in applying visibility 
constraints. 

97 At step 628, the shock scaffold engine 110 determines whether all the clusters 
have been paired. If not, the shock scaffold engine 1 10 returns to step 624, and selects 
another cluster. In another alternative embodiment of the invention, the shock scaffold 
engine 1 10 searches for a pair for the selected cluster in a layered spiral around the 
selected cluster. The shock scaffold engine 1 10 starts searching for cluster pairs in the 
closest surrounding clusters to the selected cluster. The shock scaffold engine 110 then 
visits iteratively other clusters in layers which are increasingly farther away from the 
selected cluster. Layers are organized as "onion peels" which the shock scaffold engine 
110 visits iteratively in a spiraling manner. As the shock scaffold engine examines 
clusters farther and farther away from the selected cluster, the chances of finding a pair to 
the selected cluster becomes less and less likely. When the shock scaffold engine 1 10 
finds a cluster that is invisible from the selected cluster, it eliminates new (unvisited) 
clusters associated with the invisible cluster from the next search layer. This means that 
the successive layers after the first layer become more fragmented with each iteration. 
The set of feasible cluster paths from the selected cluster therefore take the form of a 
multi-dimensional star shape (2D or 3D) having "arms" that become finer as they grow 
longer with each iteration layer. Eventually, the arms stop growing because no more 
clusters are visible or the clusters are exhausted. 
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98 At step 628, if all the clusters have been paired, the shock scaffold engine 110 
goes to the validation process, step 630, described below with regard to Figure 1 1 . 

99 In the alternative embodiment of the invention where the initial set of data 
elements is incomplete, new data elements are added over the course of deriving the 
shock scaffold, step 612 of Figure 9. At step 632, the shock scaffold engine 1 10 locates 
the new data elements in the defined clusters and proceeds to step 622 where the shock 
scaffold engine 1 10 pairs the input elements within each cluster. Such new generators 
may invalidate previously computed shocks (i.e., make their associated contact sphere 
non-empty). This effect, however, is only local and can invalidate only shocks of 
generators which will be paired with the new generators. 

100 Figure 1 1 is a flow chart of a procedure to validate shock candidates performed 
by the shock scaffold engine in accordance with the Lagrangian computational method, 
step 604 of Figure 9 and continuing from step 630 of Figure 10. The generator pairing 
process described above with regard to Figure 10 generates shock candidates, that is, 
shocks that are not necessarily valid. The shock scaffold engine 1 10 validates the shock 
candidates in order to derive the shock scaffold. A valid shock candidate has no other 
generators in its contact sphere other than the pair of generators giving rise to the shock 
candidate. 

101 At step 660, the shock scaffold engine 110 selects a shock candidate to validate. 

102 At step 662, the shock scaffold engine 110 determines whether the contact sphere 
of the shock candidate is entirely contained within the defined cluster containing the 
shock candidate. The defined cluster is one of the clusters defined at step 620 (Figure 
10). 

103 At step 664, if the contact sphere is entirely within the defined cluster, the shock 
scaffold engine 1 10 determines whether the contact sphere contains any generators not in 
the generator pair giving rise to the selected shock candidate. 

104 At step 666, if the shock scaffold engine 1 10 finds generators other than the 
generator pair in the contact sphere, the shock scaffold engine 110 defines the selected 
shock candidate as invalid. 
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105 At step 668, if the shock scaffold engine 110 does not find generators other than 
the generator pair, then the selected shock candidate is valid and the shock scaffold 
engine 110 adds it to the shock scaffold as node. 

106 At step 670, if the contact sphere is included in clusters in addition to the defined 
cluster containing the shock candidate, the shock scaffold engine 110 examines each 
cluster containing a portion of the contact sphere for generators. 

107 At step 672, the shock scaffold engine 1 10 determines whether the contact sphere 
contains any generators not in the generator pair giving rise to the selected shock 
candidate. 

108 At step 666, if the shock scaffold engine 110 finds generators other than the 
generator pair in the contact sphere, the shock scaffold engine 110 defines the selected 
shock candidate as invalid. 

109 At step 668, if the shock scaffold engine 110 does not find generators other than 
the generator pair, then the selected shock candidate is valid and the shock scaffold 
engine 1 10 adds it to the shock scaffold as node. 

1 10 After validation of shock candidates is complete, the shock scaffold engine 110 
returns to step 606 to complete the process of deriving the shock scaffold described above 
with regard to Figure 9. 

111 In the Eulerian computational method for deriving a shock scaffold, the shock 
scaffold engine 110 imposes a reference grid on the input elements and simulates wave 
propagation on the input elements through the reference grid. The reference grid 
tessellates the space containing the input elements. The logical wavefront is tracked to 
find a minimum distance to a generator forming a pair with the wavefront initialization 
generator. In a first embodiment of the Eulerian computational method, the computation 
is performed on input from a two dimensional object and the grid defines squares. In an 
alternative embodiment of the invention, the computation is performed on input from a 
three dimensional object and the grid defines cubes. In further alternative embodiments, 
other shapes are defined by the grid. Any shapes, including irregular shapes, that join 
without forming gaps (i.e. that tessellate) may be used for the grid. It should be 
understood that the method can also be equivalently executed in three dimensions. 
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Waves are simulated by sending beams from an initial grid cell in two dimensions or 
from an initial chamber in three dimensions. 

1 12 Figure 14 illustrates the three types of beams. Figures 14A, 14B, and 14C each 
show a sphere having a generator at the center 802. Figure 14A shows the first type of 
beam, a one-dimensional (ID) beam 804. The ID beam 804 propagates along grid 
directions. Figure 14B shows the second type of beam, the two-dimensional (2D) beam 
806 that propagates as a plane along two grid directions. Figure 14C shows the third type 
of beam, the three-dimensional (3D) beam 812 that is formed by combining 2D beams 
808, 810 by pair covering orthogonal directions. The ID beams 804 use paths directly 
supported by the grid and avoid overlap within the Eulerian grid. ID beams 804 however 
do not reach all the cells. 2D beams 806 fill in some of the gaps left by the ID beams 
804 in their planes of propagation. The 3D beams 812 fill in gaps left between the 2D 
planes of propagation. These beams 804, 806, 812 are referred to as Eulerian automata 
because they evolve along the grid. 

1 13 The shock scaffold engine 110 uses generators as sources of beam propagation 
and where the beams collide, shocks are formed. Each successive cell along the beam 
path is labeled with the label of the initiating generator and the signed distance vector 
from the generator. Essentially, the shock scaffold engine 1 10 derives a shock scaffold 
from input elements by initiating all types of beams from all possible generators in the 
input data set. 

1 14 Figure 1 5 is a two dimensional grid showing beam propagation using the Eulerian 
computational method of deriving a shock scaffold. Beam propagation is shown in two 
dimensions for simplicity. The 3D beam is similarly propagated through the grid. The 
two dimensional reference grid 820 has a plurality of cells, one cell 822 having three 
generators 824. Vertical ID beams 826 and horizontal ID beams 828 propagate from 
each corner of the cell 822. A 2D beam 830 propagates from one corner of the cell in a 
direction indicated by the arrow 832. Where the beam 830 passes through a node 834 of 
another cell, more ID beams 836, 838 are initialized. 

115 The shock scaffold generator 110 propagates logical waves through the grid by 
initializing beams 826, 828, and 830 at the cell 822 containing at least one generator 824. 
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The logical waves propagate through the fixed grid 820 initializing more beams 836, 838 
that carry generator and distance information in a vector. Eventually, propagating beams 
in the fixed grid collide forming a shock. 

1 16 Figure 16 is a flow chart of a procedure that is performed by the shock scaffold 
engine 1 10 (Figure 1) when operating in accordance with the Eulerian computational 
method of generating a shock scaffold. 

117 At step 900, the shock scaffold engine 110 receives as input unorganized input 
elements such as an unorganized point cloud as described above for the first 
computational method. 

118 At step 902, the shock scaffold engine 110 establishes a reference grid around the 
input elements. For a two-dimensional (2D) object, the shock scaffold engine defines a 
two-dimensional reference grid. For a three-dimensional (3D) object, the shock scaffold 
engine defines a three-dimensional reference grid. Each cell (in 2D) or each chamber (in 
3D) of the grid contains zero or more input elements, or generators. 

119 At step 904, the shock scaffold engine 110 selects a cell (or a chamber) and 
initiates vectors (that is ID beams) along the grid in a subset of possible directions of 
flow from the input elements contained in the cell. The vectors define cellular automata. 
Each vector on the grid is a cellular automaton. A cellular automaton is an entity in an n- 
dimensional grid that carries information. Each cell in the grid has a discrete state. The 
cellular automata follow rules that describe the state of the cell for the next step 
depending on the states of the cells in the neighborhood of the cell. In the present 
invention, cellular automata carry information about its generator and a distance from the 
generator. Each input from a cell along a propagation path carries information away 
from the initial generator and into the grid. ID cellular automata (that is ID beams) 
move but typically are not able to cover all grid points. The shock scaffold generator 110 
combines ID beams into 2D beams and also ID and 2D beams into 3D beams (also 
referred to as cones) that are able to propagate fully through the grid. The shock scaffold 
engine 110 continues to select cells and initiate beams until all the generators in the input 
data have been used to initiate beams. 
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120 At step 906, the shock scaffold engine 110 detects shocks wherever beams (or 
cellular automata) collide in the grid. Because beams are initiated at all generators in the 
input data, beams are propagated in many different directions in the grid and therefore at 
least some of the beams are bound to collide. A collision forms a shock which is part of 
the shock scaffold being derived from the unorganized input elements by the shock 
scaffold engine 110. 

121 At step 908, the shock scaffold engine 110 uses a Lagrangian method for finding 
shock sheets and shock curves of the shock scaffold. 

122 To trace shock waves for planar sheets and curves once the shock scaffold engine 
1 10 detects shock sources there are two additional classes of automata to propagate 
through the grid. The two additional classes of cellular automata are referred to as 
Lagrangian because they evolve away from the reference grid. The first additional class 
of automata is a 2D class of automata moving along intercepts with chambers' edges and 
the second class of automata is a ID class of automata moving along straight lines, 
keeping record of their intercepts with chambers' faces. 

123 Lagrangian cellular automaton for shock sheets, which, once initiated from a 
shock source, grows a planar sheet radially around that shock source. This 2D growth is 
obtained by computing intercepts of the sheet with chambers' edges. A ID Lagrangian 
cellular automaton for shock curves, which, once initiated from an shock source, grows a 
straight curve in one or two opposite directions from that shock source. The ID growth is 
obtained by computing intercepts of the curve with chambers' faces . 

124 Figure 17A is a diagram of a grid chamber 920 illustrating passage of the first 
type of Lagrangian cellular automata, the 2D Lagrangian automata for shock sheets. The 
grid chamber 920 is one chamber in a larger reference grid. In this example, a shock 
sheet 926 enters the grid chamber 920 from the bottom at a pair of edge intercepts 922, 
propagates through the grid chamber 920 and intercepts at the top via two more edge 
intercepts 924. Inward flows to the grid chamber 920 are indicated via arrows 926 at the 
edge intercepts 922. Outward flows to surrounding chambers are indicated via arrows 
928 at the edge intercepts 922, 924. The outward flows are used to propagate further the 
Lagrangian cellular automata for shock sheets. 



Attorney Docket No.: BRU02-01 



-so- 
ns Figure 1 7B is a diagram of a grid chamber 950 illustrating passage of the second 
type of Lagrangian cellular automata, the ID Lagrangian automata for shock curves. The 
grid chamber 950 is one chamber in a larger reference grid. In this example, a shock 
curve 952 enters the chamber from the right at a first face intercept 954 and leaves it from 
the top at a second face intercept 956. The shock curve 952 at the first intercept 954 has 
an inward flow into the grid chamber 952 indicated by arrow 958. The shock curve 952 
at the second intercept 956 has an outward flow indicated by arrow 960. The outward 
flow propagates the ID Lagrangian automaton for the shock curve 952 to other chambers 
in the grid. 

126 Figure 18 is a flow chart of the operation of shock scaffold engine 110 using the 
Lagrangian cellular automata to recover the shock sheets and shock curves of a shock 
scaffold, during step 908 of Figure 16. In general, the shock scaffold engine 110 traces 
shock waves (or logical wavefronts) by following "particles" and their intercepts with the 
fixed grid, also referred to as the Eulerian grid. The connectivity amongst these intercepts 
is then called a Lagrangian grid, tracing the flow of the particles in space. The shock 
scaffold engine 1 10 processes the Lagrangian automata in a grid chamber after the 
Eulerian wavefronts have passed entirely through the chamber. Then, for each chamber, 
the shock scaffold engine performs the Lagrangian automata computations related to that 
chamber. 

127 At step 1000, the shock scaffold engine 1 10 detects shock sheet sources within a 
selected grid chamber. The shock scaffold engine 110 collects all generator labels 
present at the eight vertices of the chamber. The shock scaffold engine 1 10 then detects 
the shock sheets sources by determining shock candidates whose contact sphere is located 
within the chamber's limits. The shock scaffold engine validates shock sheet candidates 
for those generators whose logical wavefronts have reached the selected grid chamber. 

128 At step 1002, the shock scaffold engine 1 10 performs a local iterative shock flow 
of sheets and curves within the chamber's limits. In this step, the shock scaffold engine 
110 iteratively intersects sheets to detect new shocks curves, and iteratively intersects 
shock curves to detect new shocks until all flows in the chamber are exhausted as 
described above in step 608 and 610 of Figure 9. In essence, the shock scaffold engine 
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110 performs a mini-Lagrangian shock scaffold computation (the first computational 
method described above) at the chamber scale. The shock scaffold engine 110 performs 
the local Lagrangian shock scaffold computation for all chambers containing shocks. 

129 At step 1004, the shock scaffold engine 1 10 propagates the Lagrangian cellular 
automata to other chambers and initiate new Lagrangian cellular automata from newly 
detected shocks. The shock scaffold engine 1 10 propagates the Lagrangian cellular 
automata by computing their new intercepts at the chamber's edges and faces, and by 
validating the intercepts. 

130 At step 1006, the shock scaffold engine 1 10 initiates new beams at the eight 
chambers bounding the selected grid chamber, if necessary. The shock scaffold engine 

1 10 verifies if the new beams' intercepts can minimize distance vectors, d, associated to 
the selected grid chamber's vertices. If a new beam does reset a minimum distance vector 
(that is, the beam does not collide with another), then the shock scaffold engine 1 10 
initiates a new beam to propagate further a generator label and distance function, thereby 
producing an exact distance map. 

131 In alternative embodiments of the invention, other types of data input can be used 
to build shock scaffold representations of objects. For example, a Voronoi diagram input 
can be used. Figure 19A is an example of a Voronoi diagram 1020 for a set of eleven 
point generators 1022. Figure 19B is a medial axis 1040 for the set of point generators 
1022 of Figure 19 A. Figure 19C is the shock scaffold 1060 for the point generators 1022 
of Figure 19 A. Figures 19 A, 19B and 19C each show a plurality of sampled points, the 
set of point generators 1022, taken from an object surface. Each point 1022 is contained 
in a Voronoi region 1026 defined by Voronoi vertices 1024 (also referred to as shock 
vertices) and Voronoi edges 1028. A Voronoi diagram 1020 is a decomposition of space 
into ownerships by the points 1022. The Voronoi diagram 1020 establishes 
neighborhoods by indicating space associated with each point 1022 with the boundaries 
1024. The Voronoi diagram 1020 when used as input to the shock scaffold engine 110 
indicates which points 1022 are neighbors and therefore provides rudimentary 
information for finding shock candidates. In an alternative embodiment of the Voronoi 
diagram 1020, curve segments or surface patches are used instead of points. 



Attorney Docket No.: BRU02-01 



-32- 

132 As can be seen in Figure 19C, a Voronoi graph is included in the shock scaffold 
1060. Therefore, the shock scaffold of an object surface, for example, can be used to 
compute the Voronoi diagram of the object surface. In addition, other subsets of the 
shock scaffold can be used to compute other important graphs. For example, a set of 
shock sources for sheets obtained in a first phase of the Lagrangian method indicate all 
the pairs that can be linked by Delaunay edges. A Delaunay edge is an edge such that a 
sphere with a diameter defined by the edge linking a pair of generators is empty of other 
generators. Delaunay edges are a subset of a Delaunay mesh which is a dual to the 
boundaries of the Voronoi diagram. Thus, a shock scaffold calculation can be used to 
determine a Delaunay mesh. 

133 A Gabriel graph in three dimensions is obtained by considering a set of shock 
curve sources obtained by pairing shock sheets without the need to consider shock 
vertices. A Gabriel graph is a collection of (N-l-D) simplices showing connectivity. In a 
three-dimensional case, a simplex is a triangle. Each simplex, for example, each triangle, 
must be such that the sphere passing through the points defining the simplex is empty of 
other generators. That is, the sphere passing through the three points defining the triangle 
is empty of other generators. The sphere is associated with the sources of shock curves. 
Accordingly, the Lagrangian computational method described above can be used to find 
the sphere. Accordingly, the shock scaffold computation can be used to find a Gabriel 
graph directly whereas other method of finding the Gabriel graph, such as the Voronoi 
diagram and Delaunay tessellations are indirect methods. The Voronoi diagram and 
Delaunay tessellations both require a larger structure to be computed and then searched to 
find the Gabriel graph as a subset. 

134 It is to be understood that the above-described embodiments are simply 
illustrative of the principles of the invention. Various and other modifications and 
changes may be made by those skilled in the art which will embody the principles of the 
invention and fall within the spirit and scope thereof 



